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Abstract. We analyze the evolution of the particle spectrum of the Tricritical Ising 
Model by varying the couplings of the energy and vacancy density fields. The particle 
content changes from the spectrum of a supersymmetric theory (either of an exact or a 
spontaneously broken supersymmetric theory) to the spectrum of seven particles related 
to the underlying E-j structure. In the low temperature phase some of these excitations 
are topologically charged particles that are stable under an arbitrary variation of the 
parameters. The high and low temperature phases of the model are related by duality. 
In some regions of the two couplings there are also present false vacua and sequences of 
bound states. In order to study the non-integrable features of this model we employ the 
Form Factor Perturbation Theory and the Truncated Conformal Space Approach. 
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1. Introduction 

The two-dimensional Tricritical Ising Model (TIM), which is the second in the unitary 
series of conformal minimal models after the critical Ising model p], describes critical 
phenomena in a variety of systems with tricritical points. It also exhibits superconformal 
symmetry, being the first member of the series of superconformal minimal models [2J, and 
it has interesting symmetries related to the coset constructions su{2) 2 ® sm(2) 1 /sm(2) 3 
and e(7)i <S> e(7)i/e(7) Hj. Due to these properties the TIM has attracted interest for 
many years. 

In the present paper we study the particle spectrum of the TIM perturbed by those 
relevant scaling fields which are invariant under spin reversal. These fields are the energy 
density e and the vacancy density t, with conformal weights given by A 2 = 1/10 and 
A4 = 3/5 respectively^!]. The action of the perturbed model can formally be written as 

A = A + g 2 J e(z,z)d 2 z + gi J t(z,z)d 2 z, (1.1) 

where Aq is the action of the TIM at its conformal point, whereas g 2 and #4 are the two 
coupling constants of our theory. Their dimensions are fixed by the conformal weights of 
the conjugate fields, g. t ~ A4 2 ~ 2Ai , where M. is an arbitrary mass scale of the off-critical 
model. The ratios of the particle masses of the model (II. ip . being universal quantities, 
depend only on a single dimensionless combination 77 of the two coupling constants g 2 and 
#4, that we choose as 

V= ^kt = —r ■ (1-2) 



02 



02 



In the following we write 77+ if g 2 > and 77_ if g 2 < 0. 

The evolution of the particle spectrum by varying rj is an important characteristic 
of the class of universality associated with the TIM. It is worth stressing that the action 
(11. ip provides the simplest example of a bosonic theory having kinks in its spectrum which 
do not get confined under continuous variations of the coupling constants. To appreciate 
this feature, one has to recall that, in a generic bosonic theory, the kinks are quite fragile 
objects: they get generally confined under a small perturbation of the parameters (see, 
for instance, refs. P, [7]). The stability of the kinks can be guaranteed, for instance, by 
supersymmetry but in theories that have also fermions [9]. In purely bosonic theories, 
the robustness of the kink states may only come via a fine-tuning mechanism or special 
symmetry of the perturbing operators. In the model analyzed in this paper, the latter 
property is the spin reversal invariance of both perturbations. The study of the particle 
spectrum of this model is, in any case, interesting in itself, since the action (11.11) is non- 
integrable if both g 2 and #4 are non-zero. 

1 For further details on the model and its symmetry we refer the reader to the article [5]. 



The particle spectrum of the TIM with spin reversal symmetric perturbations 



3 



2. A quick glimpse into the model 

In the TIM, the parameter gi can be interpreted as the difference between the actual 
temperature and its critical value. The #2 = critical line divides the phase plane into 
the high temperature {gi > 0) and the low temperature (#2 < 0) halves, which are 
mapped into each other by a duality transformation under which (<72,<74) <-> (— 92,9a)- 
Along the critical line g<i = 0, the model (11 .11) is integrable and has also supersymmetry 
[TU1 [TTT [12"] : the g± > half line consists of a massless renormalization group flow to the 
Ising critical point in which the supersymmetry is spontaneously broken, while the g± < 
half line consists of first order phase transition points with unbroken supersymmetry. The 
meeting point of these half lines is the tricritical point. The model is also integrable along 
the #4 = line where its particle spectrum and S-matrix are related to the E 7 Lie-algebra 
[TBI [Hj, therefore we shall refer to this line as the ^-related line. 

It is well known that the 6 Landau-Ginzburg theory has a tricritical point which is 
described by the TIM |15j . in particular the model (II. ip admits another description in 
terms of the action 



with the identification e ~ : <p 2 :, t ~ : <f) A :. In this formulation the tricritical point 
corresponds to 02 = 0,4 = 0; the spin reversal transformation acts on as 1— > —(f). This 
correspondence allows one to associate the Landau-Ginzburg potentials with the various 
parts of the (#2, g±) phase plane, as shown in Figured! 

The shape of the potentials in Figure [JJ provides a useful guide to the nature of the 
excitations at the various points of the (g 2 , #4) plane. For instance, in the high temperature 
phase, where the spin reversal symmetry is unbroken, the potential has a unique vacuum 
and therefore the particles do not have topological charge in this phase. Moving up 
from the positive horizontal axis in the first quadrant, the curvature of the minimum 
decreases until it vanishes when we reach the positive vertical axis: the theory has a 
massless spectrum at this point. Instead, moving down from the positive horizontal axis 
in the forth quadrant, two metastable vacua start to appear and they become abruptly 
degenerate with the true vacuum once the negative vertical axis is reached. Right at this 
point the theory has kink excitations, i.e. topologically charged particles, that interpolate 
between the three vacua. 

In the low temperature phase, i.e. in the second and third quadrants, the spin reversal 
symmetry is instead spontaneously broken and the ground state is doubly degenerate, 
hence in these regions we expect that the spectrum consists of kinks and possible bound 
states thereof. Moreover, it is sufficient to look at the evolution of the potential in these 
quadrants in order to conclude that the kinks of this model are always stable excitations 
(at least those with lowest mass), no matter how we vary the coupling constants. In more 
detail, starting from the negative horizontal axis and moving up in the third quadrant, 
the barrier between the two vacua decreases until it disappears when we reach the positive 
vertical axis. Hence, we expect the masses of the kinks and their possible bound states 




(2.1) 
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to decrease along this trajectory so that they adiabatically vanish once we arrive at the 
positive vertical axis. On the other hand, moving down from the negative horizontal axis 
in the third quadrant, a metastable vacuum starts to appear in the middle of the barrier: 
this state becomes abruptly degenerate with the two existing vacua just when we reach 
the negative vertical axis. At this point, the kink excitations between the two original 
vacua break down and give rise to two different sets of topological excitations relative to 
the new threefold vacuum structure. 

In the following we are going to support the physical scenario described above by 
using a combination of analytic and numerical methods. The first method is the Form 
Factor Perturbation Theory (FFPT) [6] that allows us to get analytic information on the 
spectrum in the vicinity of the integrable lines. The second method is the Truncated 
Conformal Space Approach (TCSA) [18J, which is one of the best suited numerical 
approaches for extracting the spectrum of perturbed conformal field theories. The spin 
reversal and the duality symmetries of the TIM play also an important role in the 
implementation of both methods. Let's finally mention that the two integrable cases 
g 2 = and g 4 = were studied by TCSA in [5]; the present work can thus be regarded 
as an extension of those studies to the non-integrable (#2 ^ 0, ^ 0) domain. 
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Figure 1. Landau- Ginzburg potentials associated with the various parts of the phase 
space. 



3. The two methods 



In this section we briefly discuss the two techniques that have been employed to study 
the evolution of the spectrum of our model. 



The particle spectrum of the TIM with spin reversal symmetric perturbations 



5 



3.1. Form Factor Perturbation Theory 

The Form Factor Perturbation Theory permits to investigate the behaviour of a theory 
in the vicinity of an integrable direction, either if the integrable theory is massive [6, 7jj 
or massless [8]. Let *4.; nt be the action of the integrable model and T(x) the field that 
moves the system away from integrability, so that nearby the action of the theory can be 
written as 

A = Ant + A J d 2 x T{x) . (3.1) 

If the original integrable model has scalar excitations^ Ai(6), A 2 (6), ■ ■ ., the first order 
correction in A to their mass M k is given by 

5M 2 k ~ 2AF fc T (i7r) , (3.2) 

where 

Fj(6 12 ) = (0|T(0)|A fc (^)A fc (^ 2 )) (3.3) 

is the two-particle form factor of the operator T(x), with 12 = Q\ — 9 2 . A similar formula 
also holds if the integrable model has massless excitations, with a proper interpretation of 
the rapidity variable 6 [8]. If, instead, the integrable theory has topological kinks \K ab {0)) 
interpolating between the vacuum states \a) and \b), the relevant formula is 

5M 2 ab ~ 2AF a T b (ivr) , (3.4) 

where 

Fj b (6 12 ) ee fc (0|T(0) [K^) K ba (9 2 )) . (3.5) 

In integrable field theories the form factors of a generic scalar operator T(x) can be 
computed exactly [HI [17] and therefore we can study the evolution of the spectrum near 
to the integrable lines. However, one has to be careful when applying the formulas (13.31) 
and (I3.5P since the two-particle form factor can be singular at 9\ 2 = m. This happens 
if the perturbing operator T(x) has a semi- local index 7 with respect to the particle 
excitations (this is usually the case of the topological kinks, for instance). In this case 
the point #12 = ±i7r is the location of a simple pole of the two-particle form factor, whose 
residue is given by [TJ [TTj 

-iRes ei2=±i7r F a "^ 12 ) = (1 -e^) k (0\T\0) k . (3.6) 

The result of the brief analysis given above can be summarized as follows. If Fi(m) 
is finite, the corresponding particle or kink excitation survives the perturbation and its 
mass is adiabatically shifted from the original value. Vice versa, if F l (ivr) is divergent, 
the corresponding excitation no longer survives as asymptotic particle of the perturbed 
theory, i.e. it gets confined. 

2 The rapidity variable 9 expresses the dispersion relation of the excitations. If we consider a massive 
particle of mass to, its energy and momentum are given by E = mcosh#, p = TOsinh#. If, instead, the 
particle is massless, it can be a right or a left mover excitation. The dispersion relations of a right mover 
are En = p = Me 6 , while those of a left mover El = —p = Me~ e , where M is an arbitrary mass scale. 
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3.2. The truncation method 

The Truncated Conformal Space Approach (TCSA) [18] can be used to study perturbed 
conformal field theories. It consists in a finite size analysis of the spectrum of a theory 
on a cylinder of circumference R. This is realized by truncating the conformal basis at 
a certain level n in the Verma modules of the irreducible representations of the primary 
fields and solving the eigenvalue and eigenvector problem of the truncated Hamiltonian 
operator numerically. 

To extract a quantity F relative to the R — * oo limiting model, one calculates F(R, n) 
(i.e. F at a cylinder circumference R and truncation level n) for several values of R, and 
possibly also for several values of n. Usually there is a range i?F,mm < R < RF,max(n) 
such that in this range F(R, n) is an approximately constant function of R, and its value 
can be regarded as a good approximation of the infinite-volume value of F. This range is 
called the physical window for F. Both Rp >m i n and RF,max( n ) depend on F, and RF,max( n ) 
also depends on n and is expected to tend to infinity as n is increased. Hence, the 
approximation to F usually improves if n is increased. For R > RF,max(n) the truncation 
effects spoil the approximation to F and give rise to an unphysical window of values for 
this quantity. 

We determine the masses of the particles from the energy differences 

AEi(R, n) = (Ei(R, n) - E (R, n)) , (3.7) 

where Eq(R, n) is the lowest energy line. In our case, the individual energy levels Ei(R, n) 
diverge as n — > oo if g± ^ 0, since the conformal weight of the perturbation t is greater 
than 1/2, the energy differences Ei(R,n) —Eo(R,n) are however convergent [19]. We also 
use the TCSA to calculate form factors and vacuum expectation values. This can be done 
by computing numerically the eigenvectors belonging to the various energy eigenvalues. 
The Hamiltonian operator of the model (11.11) on a cylinder of circumference R is 

R \ 12 

pR pR 

+27rg 2 \ 2 R- 2A2 / e{x, 0) dx + 27r# 4 A 4 ir 2A4 / t{x,0)dx, (3.8) 
Jo Jo 

where c = 7/10 is the central charge of the TIM, whereas L and L are the zero index 

generators of the chiral Virasoro algebras. For convenience, we have chosen to introduce 

the quantities 

A 2 = (2tt) 2A2 " 1 x 0.0928344 , A 4 = (2tt) 2A4 " 1 x 0.1486960 . (3.9) 

These are precisely the values of the coupling constants that ensure a mass gap M = 1 
when we specialize the above Hamiltonian to the integrable massive cases (g 2 = ±1,<74 = 
0) and (g 2 = 0,^4 = —1) [2D] . Further data which are necessary to implement the TCSA 
for the TIM can be found in [5]. 

To study the mass spectrum it is sufficient to consider the zero momentum subspace 
only. Moreover, one can also treat the parts of the Hamiltonian which are even or odd 
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with respect to spin reversal separatelj|§ . In our calculation we used periodic boundary 
conditions, which selects the sector of zero topological charge in the Hilbert space of the 
multi-particle states. As we will show below, there may be other phenomena related to 
the boundary conditions and to the finite size of the system that have to be properly 
taken into account in order to interpret the spectrum correctly. 

In our numerical calculations we varied the coupling constants g 2 and #4 along the 
square with corners (±1,±1) shown in Figure [2j The truncation level that we used is 
n = 8, that corresponds to a Hamiltonian truncated to 1624 conformal states. The finite 
volume spectra shown in Figures HI El [9] and [10] are always plotted as functions of mR, 
where m is the mass gap, whose value depends however on the region of the coupling 
constants which each figure refers to. 



(-1,1) 



(1,1) 



(-1,-1) 



(1,-1) 



Figure 2. The square in the (52,54) plane along which the TCSA calculations were 
done. 



4. The spectrum in the high temperature phase 

In this section we describe our results on the particle spectrum of the model (II -ip in the 
high temperature phase, made up of the first and the fourth quadrant of the (g 2 , g^j plane. 
This phase is described by the variable 77+ G (—00, +00): the value 77+ = —00 corresponds 
to the negative vertical axis, 77+ = to the positive horizontal axis, and 77+ = +00 to the 
positive vertical axis. 

In this phase the Z 2 spin reversal symmetry is exact and there is a unique vacuum 
of the theory (although in the fourth quadrant there are also two false vacua). Hence 
here the excitations are ordinary scalar particles. In order to follow the evolution of the 
spectrum it is convenient to start our analysis from the vicinity of the negative vertical 
axis, i.e. from the fourth quadrant. 

3 In addition to the two spin reversal even fields e and t, the TIM possesses other two even fields, 
the identity operator / and the supersymmetry generator G, with conformal weights and 3/2 
respectively. The spin reversal odd sector is generated by the magnetization operator a and the sub- 
leading magnetization a', with A a = 3/80 and A CT < = 7/16. 
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Figure 4. The lowest nine energy levels Ei, i = 0...8 in the first-order transition point 
as functions of mR, obtained by TCSA. Red lines are even, blue (dashed) lines are odd 
under spin reversal, m denotes the mass gap. 

4-1. The fourth quadrant 

Precisely at 7] + = — oo, the model has a first order phase transition. This means that 
along the negative vertical axis the theory has three degenerate ground states |0), |+) 
and |— ): the vacuum |0) is even under the spin reversal symmetry Q, while |+) and 
|— ) are mapped into each other: |— ) = Q\+) (see Figure Ela). Right at r/ + = — oo, 
the theory presents also an exact supersymmetry: the particle spectrum consists of four 
kinks \K_q), \K _), \K 0+ ) and \K +0 ) of equal mass. These kinks, that do not have bound 
states, provide an irreducible representation of supersymmetry and their exact S-matrix 
was obtained in [11]. The action of the spin reversal symmetry Q on the kinks is given 
by Q\K- Q ) = \K +0 ), Q\K ^) = \K 0+ ). 

On a cylinder of circumference R, the three degenerate vacua give rise to three 
exponentially split energy lines, for the phenomenon of tunneling that occurs at a 
finite volume. With periodic boundary conditions, there are however no energy lines 
corresponding to the one-particle states of the kinks. Above the three exponentially 
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degenerate lines of the ground states, one finds instead doublets of lines corresponding to 
the neutral kink-antikink configurations (see Figure H]). The proper counting of the states 
appearing in the finite volume with periodic boundary conditions is explained in Figure 
Ob. 

It is now important to understand what happens if we deform the theory away from 
this integrable situation, moving into the fourth quadrant by means of the perturbation 
Q2 J e(x)dx (g 2 > 0). The drastic effect of this perturbation is to lift the degeneracy of 
the three vacua: in the perturbed theory, the central vacuum |0) becomes the true ground 
state of the theory, whereas the other two vacua |±) become metastable ground states, 
separated from |0) by a gap 

AE ~ g 2 ((ft\e\Q) — (0|e|0)) , (4.1) 

where \Q) — |±). As a consequence, we expect that the kinks of the unperturbed system 
will get confined as soon as we move away from 7] + = — oo by switching on the coupling g 2 
of the energy operator e. This can be directly seen by the FFPT, reported in Appendix 
C: the two-particle form factors on the kink states of the e operator have in fact a pole 
at 9 = in [23]. The linear confining potential between the constituents of the two pairs of 
neutral kink-antikink states \K _K_ ) and \K 0+ K +0 ) gives rise then to a dense sequence 
of bound states, with the number of bound states going to infinity as g 2 — > [7j (see 
Figure [5]) . 



X 



A E 




(a) (b) 

Figure 5. (a) Kink-antikink state at distance x, where AE is the gap of the unbalanced 
vacua. The kink and the antikink are subject to a linear confining potential (b), with 
slope given by AE, and they give rise to an infinite sequence of bound states. 



This analytic result is indeed confirmed by our numerical calculations. In particular, 
Figure [6] shows that the number of particles below the stability threshold increases to 
infinity as the first order transition line is approached (i.e. r] + —>■ — oo). This scenario is 
similar to the one observed in the scaling region of the critical Ising model near the thermal 
axis once it is perturbed by a magnetic field [6] and it realizes the so-called McCoy- Wu 
scenario [22J. An important difference from the Ising model, however, is the Z 2 symmetry 
of our model. At r/ + = —00, the two-kink states \K _K_ ) and \K 0+ K +0 ) (the rapidities 
of the kinks are suppressed) are mapped into each other by spin reversal, so they form 
degenerate pairs (although in finite volume this degeneracy is lifted). The two-kink states 
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\K_ Kq_) and \K +Q K Q+ ) also form similar degenerate pairs, and so do the two-kink states 
\K_ K 0+ ) and \K +0 K ^). The degeneracy is lifted as soon as the perturbation is switched 
on; the arising bound states are singlets with definite parity. As the kinks do not have 
bound states at r/ + = — oo, all the bound states in the neighbourhood g 2 > of the 
transition point arise from the free two-kink configurations mentioned above. Since the 
|±) vacua are lifted, the original two-particle states 1^+0^0+); \K-qKq + ), and 

\K +0 K _) of the continuum disappear entirely from the spectrum. In their place, there 
is a sequence of bound states \B„), which arises from the even and odd superpositions of 
the \K _K_ ) and \K 0+ K +0 ) states: 

± „ | A -o-/f-o)± |Ao+AW 
y/2 

The superscript in \B^) denotes parity, the subscript numbers the elements of the tower of 
bound states. As we mentioned, \B+) and \B~) are not degenerate. Due to the exactness 
of the Z 2 symmetry the vacuum is expected to be even, and then the first state above 
it, i.e. the lowest one-particle state, is expected to be odd, and generally the parity is 
expected to alternate in the sequence of one-particle states. 

The numerical data that we obtained by varying r] + indeed confirm that in the high 
temperature phase the spin reversal symmetry is unbroken, there is a unique vacuum and 
the excitations are singlet scalar particles with alternating parity. The disappearance of 
the two-kink states \K_ K _), \K +0 K 0+ ), \K_ K 0+ ) and \K +0 K _) does not have a drastic 
effect on the numerical spectrum on the cylinder with periodic boundary conditions, since 
the \K_ K 0+ ), and \K +Q K Q _) are filtered out and |i^_o-F^o- ) and \K +0 K 0+ ) are identified 
with \K _K_ ) and \K 0+ K +0 ) under these boundary conditions: what happens is that, 
energy lines that were previously asymptotic degenerate, now become separated by tiny 
gaps. 

The doubly degenerate false vacuum can be seen in the i] + < domain in the TCSA 
spectra in the form of a linearly rising double line-like pattern (see Figure [9ja). The 
slope of the pattern, which equals to the gap AE between the false and true vacua, 
increases from as r] + is increased, and finally the pattern disappears from the finite 
volume spectrum as r] + approaches 0. AE can be calculated by FFPT around the first 
order phase transition line. At = — 1 the exact result to first order in perturbation 
theory is 

AE 



92 



0.3076... = 2A 2 x (2tt) ^ a 4 A4 4 A 2 , (4.3) 



where 



A 2 = 2.668319... (4.4) 

is the exact value of the vacuum expectation values of e given in [251 12S] , and the factor 
next to 2A 2 arises from various normalization factors used in the present paper and in 
[25| [26] (see Appendix B for more details on the calculation of AE). On the other hand, 
by TCSA we obtain 
AE 



92 



0.29, (4.5) 
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Figure 6. The particle mass ratios nii/mi as functions of r)+ in the high temperature 
phase in the range —1.361 < r/ + < 0.6 7 obtained by TCSA. Even states are shown in red, 
odd states in blue and with dashed line. The thresholds 1m\ and m\ + m,2 for even and 
odd particles are drawn with dashed thin black and green lines. 



which is in reasonable agreement with the exact value. 

Figure [6] shows the numerically calculated mass ratios m^/mi (mi being the mass 
of the lightest particle) as functions of rj + . Truncation errors tend to increase for large 
absolute values of rj + , therefore the figure shows only a limited range around rj + — 0. 
The masses mi and m2 of the two lightest particles are also shown in Figure El The 
stability threshold, taking into account spin reversal symmetry, is 2m\ for even particles 
and mi + m 2 for odd particles (if there are no even particles, then the threshold is 3mi). 

Hence, increasing 7] + , the number of stable excitations decreases. Around rj + = 0, i.e. 
near to the positive horizontal axis, there are four stable excitations. However, precisely 
at i] + = 0, the theory acquires three more stable particles, all of them above threshold. 
In order to understand the nature of the spectrum near to 77+ = 0, it is obviously 
more convenient to switch to the integrable formulation of the theory along the positive 
horizontal axis. 

4-2. The E-j-related line and FFPT 

As we mentioned in the introduction, the 77 = line corresponds to an integrable 
model related to the Toda field theory based on the exceptional Lie algebra E7. In 
the high temperature phase, this integrable model has a spectrum composed of seven 
scalar particles A iy i = 1 ... 7. Their masses m,, i = 1 ... 7 (see Table CD) and S-matrix 
amplitudes are known exactly [131 H3 ■ 

Note that A±, A3 and Aq are odd under spin reversal symmetry while the other 
particles are even. Only the first four particles have a mass below the stability thresholds 
2mi and mi + 777,2- The higher three particles are nevertheless stable along this axis due 
to the integrability of the theory. 




Figure 7. The particle masses mi (in red) and m 2 (in blue and with dashed line) 
obtained by TCSA along the square shown in Figure [H The section < 52 < 0.6, 
<74 = — 1 is not shown. 

Table 1. Masses and spin reversal parities of the particles on the Ej -related axis in the 
high temperature phase. 



mass 


mi/mi 


Z2 parity 


mi 




1 




m 2 = 


2mi cos(57r/18) 


1.28557 


+ 


m,3 — 


2toicos(2tt/18) 


1.87938 




7714 — 


2mi cos(7r/18) 


1.96961 


+ 


— 


4toi cos(7r/18) cos(57r/18) 


2.53208 


+ 


m 6 = 


4mi cos(4tt/18) cos(2tt/18) 


2.87938 




m-j — 


4mi cos(tt/18) cos(2tt/18) 


3.70166 


+ 



4-3. Decay processes 

Moving away from the axis r] + = 0, the expectation is that there exists a range of values 
of rj + where the four lowest particles are still stable, whereas the higher three particles 
decay in all possible channels compatible with the symmetry of the perturbation. The 
perturbation that breaks the integrability of the theory at r] + = and plays the role of 
the operator T(x) in eqn (13.11) is the field t, which is even under the Z2 spin reversal 
symmetry. This implies certain selection rules in the decay processes. So, for instance, 
the particle A 5 can only decay in the channel 

A 5 -> A x At , (4.6) 

even though the decay process A 5 — > A\ A 2 would be permitted by kinematics. Similarly, 
it is easy to check that the other possible decay processes compatible with the selection 
rule are 

A 6 -> A x A 2 
A 7 — >■ Ai At 

A 7 -> A 2 A 2 (4.7) 
A 7 -> A x A 3 
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(a) (b) 

Figure 8. Signature of the appearance of a resonance out of the integrability in the 
finite volume spectrum, (a) Integrable point, (b) Perturbation switched on. These are 
the plots of the energy levels versus mR, where m is the mass gap. 



A 7 -> A 2 A A 
A 7 -> A 1 A A 2 . 

The decay processes have a fingerprint in the finite volume spectra, namely the repulsion 
of the energy lines [6]. Observe, for instance, Figure Ela: this refers to the lowest energy 
difference lines above the two-particle threshold (in the even sector only) on the integrable 
E-j axis. One can see that there are several level crossings; in particular there is one around 
mR ~ 12 between the first and the second lines. The first line corresponds to the threshold 
state | At At), whereas the second line corresponds to the particle |A)- When we switch 
on the #4 perturbation, the crossing disappears, as it can clearly be seen in Figure |8jb. 
In particular, the repulsion of the first and second lines at mR ~ 12 corresponds to the 
decay A — > A A- 

In principle, the widths T^- of the decays A}~ — > AA/ can be determined from the 
repulsion of energy lines by a method proposed in [28J. In this approach, which was called 
Breit-Wigner method in [28], T is related to the two-particle phase shifts of the various 
two-particle levels |AA/)- In more detail, in order to measure T one can use the formula 

A6 = min8x(E) - max5 2 (E) = Ay/-/3T , (4.8) 



where 



^ dE 



(4.9) 

E=E k 

should be negative, E k is the center of the resonance and 5 (E) is defined by the integrable 
S-matrix as 



5 (E) = -i\n S(E) . (4.10) 

For the sake of clarity, we have suppressed the labels i, j, k relative to the particles. The 
symbols S\(E), 62(E) in (14.81) refer to the phase shifts of two neighbouring two-particle 
levels I A A/) and are given by 

6 a (E) = -R(E)p a (E), a = 1,2 , (4.11) 
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where R is the circumference of the cylinder and p a {E) is the momentum, defined by the 
equation 



E = y/fi + mt + y/rt + Tn} . (4.12) 

Hence, from the numerical determination of the energy levels we can measure E(R) as a 
function of R; from this measurement, using (14. lip , we can extract 6i(E) and 62(E), and 
then (j4.8p provides T. 

Using the above consideration we have studied the decay processes 

A 5 -> A 1 A 1 

at g 2 = 1 and small g±, finding the following numerical values for the decay rates 

I*~2.<£, T 6 12 ~ 0.2 -gl, (4.13) 
and therefore the universal ratio 

r? 



12 

r 5 
1 11 



0.1 . (4.14) 



This result shows that the decay of A$ is considerably slower than the decay of A 5 . 
Although the kinematical phase space factor $^f. jm . generally tends to suppress the decay 
of heavier particles in 1 + 1 dimensions [29J, in the present case this effect is not sufficient, 
however, to justify the result above: in fact, the phase space estimate of the universal 
ratio (14.141) gives in this case 



r? 2 $? 2 / m 2 (m 2 -Am 2 ) 

12 - ' ' 1 - 0.783... . (4.15) 



r fi $n y ( m 6 - ( m i - m 2 ) 2 )(ml - (m 1 + m 2 ) 2 ) 
Therefore the longer lifetime of A$ has to be attributed to the dynamics of the model. 
We plan to analyze this aspect in a future publication. 

4-4- Mass corrections 

Coming back to the stable particles of the spectrum, it is relatively easy to compute 
the correction to the masses of the two lowest particles once we move away from the E-j 
integrable axis. This can be done by combining the FFPT and the TCSA, together with 
a cross-checking between the two methods. The first order mass correction of the particle 
Ai moving away from the rj + = axis is given by 

6m t = ^-27r\ 4 F^(m), (4.16) 

mi 

where —0 2 ) = (0\t(0)\A i (9i)A i (8 2 )) is the two-particle form factor of t. As a matter 

of fact, i*ii(i7r), F2 2 (i7r), Fj 3 (i7r) and F4 4 (i7r) can be directly extracted by TCSA (see 
Appendix B and [2U [26]). At g 2 — 1 we obtained 

*w - • ^ • (417) 

km - • - -H ■ (4 - i8) 
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Similarly, from TCSA we can also compute the absolute value of the one-particle form 
factors F l 2 = (0\t(0)\A 2 (9)) and F\ = (0\t(0)\A A (9)). The result is 

t - °- 35 - °- 26 (4 19) 

27rA4 1 2tt\^ 

As shown in the Appendix A, there is an exact relation between the two-particle form 

factors Fl^iir) and F 2 t 2 {i'n') and the one-particle form factors F 2 and F\, expressed by 

F^O) = 11.397 ■ F\ + 10.3632 • F\ (4.20) 
F* 2 (ivr) = 23.1056 • F\ + 26.2404 • F{ . (4.21) 

If we assume that F\ < and F\ > 0, then (Q5j> . (OD]l and flOTD yield 

^ = ^ • ^ ^ • < 4 - 22 > 

which is in satisfactory agreement with f)4.17p . 

With the above values, we can extract the universal ratios 

Sm 2 / Snii = 0.76 , 

5m 3 /5m 1 = 0.84 , (4.23) 
Sm^/Smi = 1.35 . 

In order to get another universal ratio, let us look at the first order correction to the 
vacuum energy density e mc . By TCSA (see Appendix B) we obtain 

^ = 0.11 . (4.24) 

On the other hand, the exact value of this quantity is 



Se 9: 



4 



= 0.1130... = A 4 x (2tt) x " a 2 Ao 2 A 4 , (4.25) 

94 

which is in good agreement with the above numerical result. A4 = 3.70708... is the 
vacuum expectation value of t calculated in (251 126] (the factor next to A 4 arises from 
various normalization factors used in the present paper and in [25l 126]). This allows us 
also to extract the universal ratio 

= -0.086 . (4.26) 
4-5. The first quadrant 

After the analysis in the previous section we can now proceed to investigate the evolution 
of the spectrum in the first quadrant. First of all, note that the corrections to the two 
lowest masses have the same sign, which depends on the sign of g^. We can use the 

(2) 

above results for 5mi and 5m 2 to estimate the value of n + where the mass r«2(?7+) of the 

second particle reaches the threshold 2mi(i] + ) of the lowest particle. This happens for a 

(2) 

value of rj + > of the first quadrant, given by r] + ~ 1: although this is higher than the 
value r/ + ~ 0.7 extracted from the numerical data shown in Figure El it is nevertheless 
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a reasonable estimation of this quantity since the theoretical result was based on just 
first-order perturbation theory. 

Besides the three higher particles A 5 , A e , A 7 that were stable only at r/ + = and 
decay as soon as we move away from the horizontal axis, the TCSA analysis shows that 
the same pattern as was seen in the fourth quadrant also persists in the first quadrant: 
namely, increasing the value of r/ + , the number of stable particles continues to decrease. 
The first particle that disappears into the threshold 2mi(i] + ) is A 4 at the critical value 
77+ , which can be estimated from first-order perturbation theory to be r]± ~ 0.04. This 

is followed by A 3 that disappears into its lowest threshold, given by (mi(r] + ) +m 2 (?7 + )), at 

(3) (3) 
the critical value r] + , for which perturbation theory gives r] + ~ 0.62. These theoretical 

results for 77+ and 77+ are in reasonably good agreement with the numerical data (Figure 

[6]). An example of the finite volume spectra calculated numerically is shown in Figure 

[Tol a at r] + = 0.3, where the theory contains 3 stable particles. 

As just described, moving toward the positive vertical axis, there is a depletion of 
the number of stable particles, until there remains only one in the neighbourhood of 
i] + = +00. At the same time, notice that increasing rj + , the value of this lowest mass 
becomes also smallei0. The mass gap of the theory finally vanishes when 77+ = +00, i.e. 
when we have reached the vertical axis. 

During this evolution, however, there has also been a qualitative change of the 
spectrum: in fact, reaching the positive vertical axis the lowest excitation has turned 
into a fermion. This becomes evident by looking directly at the nature of the theory at 
7] + = +00 (notice that the positive vertical axis also corresponds to 7/_ = +00). 

4.6. The spontaneously SUSY breaking axis 

The analytic control of the theory in the vicinity of r]± = +00 is provided by the 
integrability of the model along the vertical axis, where the system is also supersymmetric. 
Along the positive vertical axis the supersymmetry is however spontaneously broken: the 
low-energy excitations are given in this case by the massless right and left mover Majorana 
fermions, which play the role of goldstino pJJ [TTJ [12] . The factorized scattering theory 
was proposed in [UJ and the basic form factors were calculated in [21]. The massless 
Majorana fermions are nothing else than those of the Ising model, connected to the TIM 
by the massless renormalization group flow that occurs along this line. 

Breaking the integrability of the r)± = 00 model by means of the operator T = e, 
the left and right moving excitations become adiabatically massive, as it can be directly 
observed in the TCSA data (see Figure [7]). To compute the mass m of the fermion 
generated by the perturbation, we need to employ in this case the massless FFPT [8] : at 
the lowest order, we have 



4 This behavior is in agreement with the first-order correction (|4.16p to mi, that is negative for 54 > 0. 



m ~ (?2 lim 




(4.27) 
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where the two-particle (right-left) form factor of the e operator is given by [21J 




(4.28) 



Using ( 14.27P and (14.281) . one can easily check that m is a finite quantity, positive for g 2 > 
and negative for g 2 < 0. Note that, for a Majorana fermion, a negative value of the mass 
signals that we are in a low temperature phase, i.e. that the theory has two degenerate 
vacua. This is the topic of the next section. 

5. The spectrum in the low temperature phase 

The low temperature phase is composed of the second and the third quadrants. This 
phase is described by the variable G (— oo, +oo): rj- = oo corresponds to the positive 
vertical axis, t]_ = to the negative horizontal axis and ?]_ = — oo to the negative vertical 
axis. 

In this phase the model presents generically two degenerate ground states |— ) and |+), 
which are mapped into each other by the spin reversal symmetry operation Q: |+) = Q\ — ) . 
An important feature that is worth stressing is the following: no matter how we vary the 
coupling constants in this parameter region, the two vacua |— ) and |+) always remain 
degenerate, the only changes being in the shape and the height of the barrier between 
them. Hence, in this model the kink states that interpolate between the two vacua are 
always stable excitations. As we will comment below, this can be explicitly checked and 
confirmed by a FFPT computation at i]_ = 0, where the theory has the E-j structure. 

The only points where one has to be careful are the two limiting cases = oo and 
?7_ = — oo. In the former case, the barrier between the two vacua disappears and, as we 
saw at the end of Section H~5l the kinks become massless. In the latter case, the two vacua 
|— ) and |+) and the false vacuum |0) that emerges between them in the region 77 < 
become degenerate precisely at r]_ = —00. Hence, at this point, the original kinks break 
into a new set of smaller kinks, i.e. those relative to the three vacua of the first order 
phase transition of the TIM. 

To obtain the spectrum of the theory in the low temperature phase we can take 
advantage of the duality of the model. At a formal level, note that the masses at a given 
point (#2? #4) in the low temperature phase are the same as the masses in the dual point 
( — (72, c/4) in the high temperature phase. This result can be obtained as follows: since e 
is odd and t is even under the duality transformation D: 



However, to interpret correctly the TCSA data, one needs to take into account the 
existence of kink states and the periodic boundary conditions imposed to the system. 



D~ L eD = -e, 
D~HD = t, 



for the Hamiltonian we have 




(5.1) 
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Note that D maps the even sector of the Hilbert space onto itself, so in the even 
sector H(g 2 ,g4) has the same spectrum as H(—g 2 ,g4). Concerning the odd sector of the 
spectrum, the particles with odd parity in the high temperature phase become instead 
kinks in the low temperature phase, i.e. topologically charged states that interpolate 
between the two vacua. However, the periodic boundary conditions that we have imposed 
on the TCSA filter only topologically neutral states. Hence, by switching g 2 — > —g 2 , 
we expect that the net result in the numerical outcomes will be the disappearing of 
the energy levels corresponding to the odd particles in the high temperature phase: in 
the low temperature phase these lines should become exponentially degenerate with the 
energy lines of the even levels. In other words, in the low temperature phase the TCSA 
data should show doubly degenerate lines and each doublet should have an even and 
odd member. This is indeed the clearly shown in Figures [9] and [TQl where we 

present the spectrum of the theory at given values (#2, #4) and (— gi-,gi)- This feature 
is nothing but a finite volume manifestation of the spontaneous breaking of the Z 2 spin 
reversal symmetry: the two vacua (and all the excitations above them) have in fact an 
exponential splitting of their energy for the tunneling phenomenon that occurs at finite 
volume. 

Summarizing the considerations above, in the low temperature phase and at infinite 
volume the model has two ground states |+), |— ), each ground state has a sequence 
\C n )± of topologically neutral particles above it, and there is also a sequence of kinks 
and antikinks |X>+_), \T> m + ) which interpolate between the vacua§. Under spin reversal 
symmetry the ground states, particles and kinks have the following property: 

l±> = Q\t) ; 

|Cn)± = Q\C T i)t i 

\V n + _) = Q\V%) . 

The masses of the neutral states at a given point of the plane of the coupling constants 
are the same as the masses of the even particles in the high temperature phase at the dual 
point, whereas the masses of the kinks are the same as the masses of the odd particles in 
the high temperature phase. In finite volume with periodic boundary conditions, however, 
we do not see individual kinks or other topologically charged configurations. We do see the 
vacua, the neutral particles and only neutral configurations of two or other even number of 
free kinks, which, moreover, get mixed in finite volume into even and odd eigenstates of the 
spin reversal operator and acquire an exponential splitting due to tunneling effects. With 
periodic boundary conditions there are also identifications between certain configurations, 
similarly as in the first order transition point (see Figure [3]). 

After the discussion above we can proceed to describe the spectrum in the second 
and third quadrants in detail. 

5 The indices n and m number the members of the sequences, while + and — refer to the vacua between 
which the kinks interpolate or which support the neutral states. 
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(b) g 2 = -1, 4 = -0.4 

Figure 9. The first 14 energy level differences Ei — Eq, i — 0...13, as functions of mR 
in the (a) J^th and (b) 3rd quadrants at the dual values of the couplings constants. Even 
levels are in red, odd levels are in blue with dashed line, m denotes the mass gap. 

5.1. The third quadrant 

At = — oo the model is at the first order phase transition point, which was described 
in Section 14.11 As we move into the third quadrant by switching on the perturbation 
(?2 / e(x)dx (02 < 0), the degeneracy of the three vacua is lifted. In contrast with the 
fourth quadrant, the central vacuum |0) now becomes a metastable ground state, and the 
two other vacua |±) remain true ground states. The gap AE separating |0) from |±) at 
a given point is equal to the gap at the dual point in the fourth quadrant. The expected 
consequence of the lifting of the degeneracy of the vacua is again the confinement of the 
kinks of the unperturbed system, which is confirmed by the same FFPT argument as for 
the fourth quadrant. The linear confining potential between the kinks gives rise to a dense 
sequence of bound states. As in the fourth quadrant, the number of these bound states 
goes to infinity as 02 — > 0. In the third quadrant, however, the two-kink states \Kq— K-q) 
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(a) g 2 = 1, g 4 = 0.3 
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(b) g 2 = -1, 94 = 0.3 

Figure 10. The first 14 energy level differences Ei — Eq, i — 0...13, as functions of mR 
in the (a) 1st and (b) 2nd quadrants at the dual values of the coupling constants. Even 
levels are in red, odd levels are in blue with dashed line, m denotes the mass gap. 

and l-Ko+^+o) are those which disappear from the spectrum, and the bound states arise 
from the other two- kink states |i^_o-^o+) ? \K+oKo-), l-^-o-^o-) and 1^+0-^0+)- Some of 
the bound states above are topologically charged, i.e. they are kinks, while others are 
neutral. The neutral particles arise from the neutral two-kink configurations: 

\C n ) + ~ |K+o^o + ), \C n )- ~ \K-oKo-). (5.2) 

It is important to note that \C n ) + and |C n )_ are degenerate for any given n, due to 
the spontaneously broken Z2 symmetry, although their degeneracy is exponentially lifted 
in finite volume. The topologically charged particles arise from topologically charged 
two-kink states: 

\V%) ~ \K_ K 0+ ), \V%) ~ \K„ K 0+ ), (5.3) 

where the subscripts — h and H — in |X>™ + ) and denote the vacua between which 

these kink configurations interpolate. |X > " + ) and \D'^_) are mapped into each other by 
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the spin reversal: Q\V™ + ) = |P? ), therefore they have equal mass. 

As we mentioned, under duality an even particle \B+) in the high temperature phase 
corresponds to the neutral particles \C n )± of the same mass in the low temperature phase, 
and an odd particle \B~) corresponds to a kink-antikink pair |X>™_), \T>™ + ) of the same 
mass. 

The neutral particles \C n )± (more precisely, their even and odd superpositions) can 
indeed be seen in the TCSA spectra, whereas the kinks and |T>™ + ) are filtered out 

by the periodic boundary conditions. The neutral two- kink states constituted by [D+J) 
and II}",) can nevertheless be observed in the numerical data. The disappearance of 
the original two-kink states |i^o— K"-o) an d \Kq + K + q) does not have drastic effect on the 
spectra with periodic boundary conditions for similar reasons as for the third quadrant. 

The false vacuum with even parity can also be seen in the ?7_ < domain in the 
TCSA spectra in the form of a linearly rising line-like pattern (see Figure [9lb). The 
duality relates this metastable vacuum to the metastable vacua in the first quadrant. 

5.2. The Ej -related line 

The number of stable particles at any given point in the low temperature phase is equal 
to the number of particles in the dual point in the high temperature phase, and this 
number decreases as r?_ is increased. The number of stable particles around = is 4, 
however at r]_ = 0, i.e. on the negative horizontal axis, 3 more particles become stable 
above threshold. These particles are related by duality to the three highest particles 
existing at rj + = 0. In particular, there are two neutral particles \C 5 )±, \C 7 )± of masses 
m 5 and m 7 (see Tabled]) which correspond to the even particles \A 5 ), \A 7 ), and there is 
one kink-antikink pair |Z?1_), \T>": + ) of mass m 6 corresponding to the odd particle \Aq). 

It was stressed in the introduction that an important property of the model is 
the presence of kinks in the low temperature phase which do not get confined under 
perturbations. At rj^ = this non- confinement can be explicitly verified by FFPT. For 
this, one has to show that the form factors giving the first-order corrections to the masses 
of \T> l + _) and \D l _ + ) are finite. This is indeed the case, since these form factors 

H*(o)iPi + (i7r)p;_(o)> , 
(+|t(o)|p;_(i7r)Pi + (o)) , 

(where % = 1,3) are equal by duality to the form factors (0\t(0)\Ai(m)Ai(0)) , and the 
latter are finite since t does not have non-trivial semi- local index with respect to the \Ai) 
particles, which are local excitations. 

The particles \C 5 )±, \C 7 )±, |^-+) are unstable near the ^_ = axis, the 

allowed decay processes can be obtained from (14. 6[) and (14.71) by replacing Ai by the 
appropriate dual particles. For example, (14.61) corresponds to 

C% -> V\_V\ + , Cl -> V l _ + V\_ . (5.4) 
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5.3. The second quadrant 

Similarly to the third quadrant, the particle spectrum in the second quadrant is related 
by duality to the spectrum in the first quadrant. Even particles in the first quadrant 
correspond to neutral particles, odd particles correspond to kinks. The number of particles 
decreases from 4 to 1 as increases from to oo, the critical values where the particles 
and kinks cross the threshold are the same as in the first quadrant, i.e. r/j 4 , r]+\ 77+ . 

(2) 

For rj- > rj + there is only one kink-antikink pair in the spectrum. Similarly to the first 
quadrant, metastable vacua are not present. As ?]_ goes to infinity, the potential barrier 
between the two vacua |±) and the mass of the kink-antikink pair decrease until they 
finally vanish at the second-order transition point r/_ = oo. An example of the finite 
volume spectra calculated numerically is shown in Figure [TUlb at i]^ = 0.3, which is dual 
to the point r] + = 0.3 that the example Figure [Tola for the first quadrant shows. 

6. Conclusions 

In this paper we have studied the particle spectrum of the TIM with spin reversal 
symmetric perturbations. This is the simplest bosonic non-integrable field theory where 
kink excitations do not get confined by changing the coupling constants over a wide range 
of values. By varying the couplings, the model interpolates between a SUSY theory (either 
in its exact or in its spontaneouly broken phase) and a theory ruled by the exceptional 
algebra E 7 , with a spectrum given by purely scalar particles (in the high temperature 
phase) and kinks and bound states thereof (in the low temperature phase). 

Our results, obtained both by the FFPT and TCSA methods, are also in agreement 
with the Landau- Ginzburg picture of the TIM. We have found that the particle masses 
are the same at dual points in the high and low temperature phases, whereas the spin 
reversal symmetry properties and the topological charges of the particles are different in 
the two phases: even particles in the high temperature phase correspond to topologically 
neutral particles in the low temperature phase, odd particles correspond to (topologically 
charged) kink states. The number of stable particles tends to infinity in the vicinity of 
the first-order transition line, which is similar to the accumulation of particles found in 
the low temperature Ising model perturbed by a magnetic field [22]. However, in contrast 
with the Ising model, in TIM the parity is a good quantum number and the corresponding 
particles emerging from the kink-antikink threshold carry even or odd parity eigenvalues. 

We have also computed the first-order corrections to the lowest four masses and the 
vacuum energy density at the ^-related line, as well as the corrections to the energy 
densities of the three degenerate vacua at the first-order phase transition line. All these 
analytic results are in good agreement with the numerical estimates extracted by means 
of the TCSA. 

Finally, we would like to mention that the three-frequency sine-Gordon model also has 
a tri critical point, the neighbourhood of which was studied by the TCSA in [30], [31] . In the 
light of the results of this paper, it would be interesting to investigate more thoroughly 
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the first order transition in the three-frequency sine-Gordon model and, in particular, 
to study the evolution of the particle spectrum of this model. Moreover, the strategy 
adopted in the present paper can be also used to study the spectrum of higher multi- 
critical theories, described by higher minimal models of conformal field theories perturbed 
by several operators. 
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Appendix A 

Form Factor Bootstrap 

In this appendix we explain the derivation of the relations (14.201) and (14.211) by the form 
factor approach. The two-particle form factor F^(9) can be written as [Tj 

Ftm = QlM-j^T> (A - 1} 

where F™ m (0) and Du(9) can be obtained by specializing the general formulas [7J to the 
present case: 

D n (9) =P 10 (6)P 2 (6) 

D 22 (9) = P 12 (0)P 8 (9)P 2 (0) (A.2) 
F™ n (9) = -ismh^g w (9)g 2 (9) 



where 



F™ n (9) = -isinh^ 12 (% 8 (% 2 (0) 



cos(Y|7r) — cosh 9 



g n (9) = exp 



00 dxcosh[(n/18- l/2)x] . 



2 / T—^sm^im - e)x/2ir] 

x cosh x/2 sinh x 



(A.4) 



The Qli(9) are polynomials in cosh(6'). An upper bound on the degree d P of these 
polynomials can be obtained from the following general formula for the asymptotic 
behaviour of form factors [7]: 

lim (0i, 9 n ) ~ e^'^', y v < A„. (A.5) 

\Oi\— >oo 

This relationship fixes dp — 1 for both and <3 22 (#) : 

Q*W = 4 + ?4cosh(0), i = l,2. (A.6) 

These coefficients a*j and 6^ can be expressed in terms of the three-particle coupling 
constants T^, rf 1( r| 2 , r 22 an d the one-particle form factors F$, F\ by means of the 
residue equations 



— i lim ( 

ab 


9 - Kb 


)FI(0) - 


= r c ,F^ 

at) c 


(A.7) 


for the fusions 






10 

71", 

18 ' 




A t x Ai - 


>A 2 , 


u 2 n = 


(A.8) 


Ax x Ax - 


> -44, 


«ii = 


2 

71", 

18 ' 


(A.9) 


A 2 x A 2 - 


>A 2 , 


M 22 = 


12 

18 ' 


(A.10) 
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g 

A 2 x A 2 -> A 4 , *4 = — 7T, (A.ll) 

where it£ b denotes the fusion angle. The expressions we found are the following: 

a\ x = 1.14107 ■ V 2 n Fl - 1.77654 • T\ x Fi (A. 12) 

6* u = -1.21431 • T^F* - 10.2307 • r^Fj (A.13) 

<4 = 0.306459 • r^Fj - 1.37383 • r| 2 F 4 * (A.14) 

6* 22 = -1.76483 ■ V\ 2 Fl - 2.74766 ■ r^F*. (A. 15) 

In terms of these coefficients 

i^(i7r)=4-i4, i = l,2. (A.16) 

The squares of three-particle coupling constants can be obtained from the S-matrix 
elements 

Sniff) = -fio(0)f 2 (9) (A.17) 

S 22 {9) = f 12 (9)fs(e)M6), (A.18) 



where 

, tanhi(6» + i7r^) 
M9) = t^fe-hr-) ' (A ' 19) 
by using the residue equation 

-ilim(^-i< 6 )S40) = (r:j 2 (A.20) 

for the fusions flA.8l) - flA.lll) . We obtained 

(r^) 2 = 4.838705173 2 (T^) 2 = 1.225805260 2 (A.21) 

(r 2 2 ) 2 = 11.15518618 2 (Y\ 2 f = 19.10015279 2 . (A.22) 

Assuming > 0, r| > 0, these results can be substituted into flA.12p - flA.15p . and then 
f lA~T6l takes the form of fl4T20]) and pill . 

Appendix B 

Calculation of form factors by TCSA 

A one-particle form factor of a field <p of conformal weights (A, A) at infinite R can be 
obtained as the limit of the finite R form factor: 

F% = (0|p(0)|A) = lim R- 2A y/R^(0\<p(0)\A) R , (B.l) 

where the subscript R on the right hand side indicates that the matrix element should be 
calculated in the theory defined at size R and m is the mass of the particle \A). In this 
formula \A) is a zero-momentum state and it is assumed that (0|0)# = (A\A)r = 1. In 
the framework of the TCSA the absolute value of the matrix elements (0|</?(0)|A)# can 
be calculated. For a general matrix element the formula 

(A\<p(0)\B) Ea, b KB b (a\<p(0)\b) 



(B.2) 
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holds, where the subscript R is suppressed, it is not assumed that the eigenvectors \A) 
and \B) are normalized to unity, A a and Bj, are expansion coefficients with respect to the 
conformal basis used in the TCSA: \A) = A a |a) , \B) = J2b^b\b), and G a b = (a\b) is 
the inner product matrix of the conformal basis vectors. The expansion coefficients A a 
and Bb are calculated numerically up to overall constant factors; they are provided by the 
routine that one uses for the diagonalization of the Hamiltonian operator. 

A two-particle form factor F AA (6) at 9 = in at infinite R can be obtained as the 
following limit of finite volume matrix elements: 

F* A (i-K) = (A(0)\<p(p)\A(0)) 

= lim R- 2 *Rm((A\<p(0)\A) R - (0^(0)|0) fl ). (B.3) 

it— »oo 

It is assumed that A is a self-conjugate particle, and normalization conditions similar to 
those for (IB. II) apply. In the framework of the TCSA the matrix elements (y4|c^(0) |^4)i? 
and (0|(^(0)|0)ij can be calculated in a similar way as the matrix element in ( IB. II) . i.e. by 
the formula (lB~2l) . 



Appendix C 

Calculation of corrections to vacuum energy densities 

The first order correction (14.241) to the vacuum energy density e vac = lim^oo E vac (R)/R 
at g 2 = 1; 94 — on the ^-related axis was calculated by the following formula: 

Se vac = 94^(94 = 0) = lim 2 m A 4J R~ 2A4 (0|t(0)|0) K . (C.l) 
dg 4 i?-+oo 

The numerical calculation of the first order corrections to the energy densities of 
the |+), |— ), |0) vacua existing at the first order phase transition point is slightly more 
complicated, since these vacua are degenerate in infinite volume (R = oo). Because of this 
degeneracy, one has to use, in general, degenerate perturbation theory, i.e. one obtains 
the first-order corrections by calculating the eigenvalues of the 3 by 3 matrix constituted 
by the 9 matrix elements (+, — , 0|if p |+, — , 0) of the perturbing operator H p between the 
three vacua. Those combinations of the vacua which evolve into the split vacua existing 
at nonzero perturbation are given by the eigenvectors of this matrix. 

More precisely, in our case the first order corrections to the energy densities of the 
vacua are given by the eigenvalues of the matrix 

Mij = lim 2irg 2 \ 2 R- 2A z(U l \e(0)\U J }R, (C.2) 

R— >oo 

where i,j = 1,2,3, 94 has the value —1, and \Ux)r, \U2)r, \Us)r are the three vacua at 
cylinder circumference R, normalized to 1. The numbers {Ui\e(0)\Uj) r can be calculated 
as described above (see equation (|B.2[) ) . It is important, in general, that \U\)r, \U 2 )r, 
I £73) # should be orthogonal, otherwise the inverse of their inner product matrix has to be 
included in (1C.2I) . In our case the vacua are split at finite R by tunneling effects, so \Ux)r, 
\U 2 )r, \Us)r belong to different eigenvalues, therefore they can be identified uniquely and 
they are orthogonal. 
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The problem of the diagonalization of (Ui\e(0)\Uj) r can be simplified by symmetry 
considerations in our case: two of the three vacua \Ui)r, \U 2 )r, \U 3 )r lie in the spin 
reversal even subspace, the third one lies in the odd subspace. We choose \Ui)r and 
\U 2 )r to be the even vacua and \U 3 )r to be the odd one. Due to the spin reversal 
symmetry of the e perturbation, the matrix elements (Ui t2 \e(p)\Us)R are all 0, therefore 
the diagonalization problem can be treated separately in the even and odd sectors. The 
odd vacuum (which becomes in infinite volume) thus evolves into the odd (possibly 

metastable) vacuum of the perturbed theory, and the first-order correction to its energy 
density is given by 

8e U3 = lim 27rg 2 \ 2 R- 2A *(U 3 \s(0)\U 3 } R . (C.3) 

R^oo 

In the even sector, one has to consider the 2 by 2 matrix (Ui i2 \e(0)\Ui i2 ) r. Since the duality 
D maps the even sector into itself and commutes with the unperturbed Hamiltonian 
operator in this sector, and \U\)r and \U 2 )r have different energies, \U\)r and \U 2 )r have 
to be eigenstates of D. The eigenvalues can be +1 or —1, because D 2 = 1. We also 
know that e changes sign under duality, therefore {U\\e{Q)\Ui) r = (U 2 \e(0)\U 2 ) r = 0. 
e is also self-adjoint, so (Ui\e(0)\U 2 ) r = (U 2 \e(0)\Ui)* R . This number can be nonzero if 
the eigenvalues of \U\)r and \U 2 )r with respect to D have opposite sign. We denote the 
eigenvectors of the matrix (Ui t2 \e(0)\Ui t2 ) r by \V+)r and |VL)_r. The eigenvalues of the 
matrix (U 1)2 \e(0)\U 1)2 )r are ±\(Ui\e(0)\U 2 ) r\, so the correction to the energy densities of 
the two vacua |V+)oo and \V-)oo are 

5e v+ = lim 27n ?2 A 2 iT 2A2 |(f/ 1 K0)|£/ 2 ) R |, (C.4) 

and 

5e v _ = lim ^ng^R-^KUMO^R}. (C.5) 

Here we have chosen |V^_)^ to have the positive eigenvalue and | VL)^ to have the negative 
one. Note that these formulas would not be meaningful without taking absolute values 
on the right hand side, since \Ui)r and \U 2 )r are defined up to arbitrary phase factors. 

The analytic result ( 1C.9I) described below, the Landau- Ginzburg picture and our 
numerical calculations of energy spectra all show that one vacuum gets a negative 
correction and two vacua get the same positive correction as we move into the high 
temperature phase. It follows then that 

5e Vi = 5e v+ = -5e v _, (C.6) 

and AE = 25eu 3 - The direct TCSA calculation of 5ejj 3 and 5ey + confirms the equality 
5ejj z = Sey + and yields the result (14. 5 p for AE. 

Having discussed the calculation of the energy density corrections, we turn to the 
description of the relation between the vacua |£/i)oo, ^2)00, ^3)00; |^+)oo, l^-^oo and |+), 
|— ), |0). Equation flC.6j) implies that 

(UslelUs)^ = |([/i|£|[/ 2 }oo|, (C.7) 
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and as we mentioned above, we also have 

(Ux\e\U 2 ) oo = {UMUi)*^ (C.8) 

and all other matrix elements of e between |Z7i)oo, 1^/2)005 ^3)00 are zero. (t/i|e|Z/2)oo can 
be made real and positive by a phase redefinition of lUi)^ or \U 2 )oo- 

The matrix elements of e between the R = 00 eigenvectors |+), |— ), |0) were 
calculated exactly in [25j 12S] : 

(+|e|+> = = — <0|e|0> = A 2 > 0, (C.9) 

and the nondiagonal matrix elements are (the numerical value of A 2 is given by (14.41) . 
but it is irrelevant here). This result, together with the spin reversal properties, (10.70 
and (1C.8I) . allows us to identify |£/i)oo, ^2)00, 1^3)00 in terms of |+), |— ), |0) as follows: 

\U 1 ) OQ = a l+)+ 2 H (CIO) 

l^oo = a ' +) t h) -^7f ( C - U ) 
\U 3 ) OQ = 1 l+) ^~\ (C.12) 

where it is also assumed that the phases of |£/i)oo and 1^2)00 are chosen so that (Ui\e\U2)oo 
is real and positive, a, (3 and 7 are unknown phase factors which are not determined 
by the information we have described so far, and they can be set to 1 by the following 
redefinition of the phases of | +), |-), |0) and \U 3 ) R : |+) -> |-) -> |0) -> i|0>, 

|^3)i? ~~ * ^|^3)r- These redefinitions leave the previously stated relations involving these 
vectors unchanged. |^+)oo and |V_)oo can now be expressed as 

l^ + )oo - ^ - — j^- (C.13) 

IVLJoo = ^ l)o ° "l 1 ^ 00 = |0). (C.14) 
v2 

We note finally that since | C^i) 00 and ^2)00 are duality eigenstates with eigenvalues of 
opposite sign, equations (IC.lOj) and (IC.lip imply that the duality maps |0) and 
into each other. 



Appendix D 
Kink form factors 

Along the first-order phase transition line (g 2 = 0, < 0), when the TIM has three 
degenerate vacua and an exact supersymmetry, in the infinite volume (R = 00) the two- 
kink form factors 

(O|e(O)|# o _(0i)#-o(0 2 )> =Ft_{e x -e 2 ) , 
(O|e(O)|# o+ (0i)#+o(02)> =F o £ + (^-0 2 ) , 
(-|e(0)|K_o(^i)^o-(^)) = Ffo(^i " 2 ) , 
(+|£(0)|K +0 (^)Ko + (^)) = F^{9 X - 6 2 ) 
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are the following 

p-5(w+M) 

FL (9) = -i(L7L - tf o )____F o (0) , (D.l) 

psmh £(0 — i7r) 

e i(7T+i0) 

Ff o (0) = i(Z7_ - Eft) — f- , F (g) , (D.2) 



where 



psinh -{9 — m) 



C7b = <0je|0> , U- = (-\e\-) = (+\e\+) 



F (9) = — isinh - exp 



dx sinh(l-p)§ sin 2 (i7r - 9)^- 
U x sinh y cosh | sinh x 



(D.3) 



7 = — ln2 , p = 4 . 
These form factors have a pole at 17T, as can be explicitly seen from (ID. II) . (ID. 21) . 
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